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ABSTRACT 

We use a suite of cosmological hydrodynamic simulations to systematically quantify 
the accretion rates of baryons into dark matter halos and the resulting baryon mass 
fractions, as a function of halo mass, redshift, and baryon type (including cold and 
hot gas). We find that the net baryonic accretion rates through the virial radius are 
sensitive to galactic outflows and explore a range of outflow parameters to illustrate 
the effects. We show that the cold gas accretion rate is in general not a simple uni- 
versal factor of the dark matter accretion rate, and that galactic winds can cause star 
formation rates to deviate significantly from the external gas accretion rates, both 
via gas ejection and re-accretion. Furthermore, galactic winds can inject enough en- 
ergy and momentum in the surrounding medium to slow down accretion altogether, 
especially in low-mass halos and at low redshift, but that the impact of outflows is 
suppressed with increasing halo mass. By resolving the accretion rates versus radius 
from the halo centers, we show how cold streams penetrate the hot atmospheres of 
massive halos at z ^ 2, but gradually disappear at lower redshift. The total baryon 
mass fraction is also strongly suppressed by outflows in low-mass halos, but is nearly 
universal in the absence of feedback in halos above the UV background suppression 
scale, corresponding to circular velocities ~ 50 km s~^. The transition halo mass, 
at which the gas mass in halos is equal for the cold and hot components, is roughly 
constant at ~ 10^^'^ and does not depend sensitively on the wind prescription. We 
provide simple fltting formulae for the cold gas accretion rate and the corresponding 
efficiency with which dark matter channels cold gas into halos in the no-wind case. 
Finally, we show that cold accretion is broadly consistent with driving the bulk of the 
highly star- forming galaxies observed at z ~ 2, but that the more intense star formers 
likely sample the high end of the accretion rate distribution, and may be addition- 
ally fueled by a combination of gas recycling, gas re-accretion, hot mode cooling, and 
mergers. 

Key words: cosmology: theory ~ galaxies: formation, evolution, high-redshift - hy- 
drodynamics 



1 INTRODUCTION 

Both observations and theoretical models indicate that the 
accretion of gas from the intergalactic medium (IGM) is a 
fundamental driver of galaxy formation and evolution. Such 
accretion is necessary to explain the observed evolution of 
cold gas in and around galaxies (e.g., Prochaska & Wolfe 
20091. Measurements of the cold gas reservoirs and of gas 



consumption time scales also imply that a continuous supply 
of gas from the IGM is required to maintain star formation 
over most of the Hubble time (e.g., Erb|2008 Bauermeister 
et al.|2010| . At the same time theoretical models of galaxy 
formation in a ACDM universe naturally predict continu- 
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ous accretion of gas into galaxies (e.g., Murali et al. 2002 
Katz et al.]|2003| [Springel fc Hernquist 2003b' IKeres et al 



et al 



2005| IKeres et al.||2009a| |Dekel et al 



20111 



2009^ i_van de Voort 



The rates at which halos accrete dark mat- 
ter have now been well quantified across cosmic time using 
A*'— body simulations, as a function of halo mass and en- 
vironment, as well as separated according to whether the 



accretion is smooth or in the form of mergers (e.g., Wech- 


sler et al.|2002 


Neistein & Dekell2008||Fakhouri & Ma|2008 


20091 


20101 Fa 


khouri et al.|2010||McBride et al.|2009| Genel 


et al. 


20101. The recent progress on this front was possi- 



ble thanks to the large dynamic range provided by A''— body 
simulations such the Millennium and Millennium II simula- 
tions ( [Springel et al.|2005[ [Boylan-Kolchin et al.|2009| ). The 
statistics measured from those dark matter simulations, and 
the derived fitting formulae, are the basis for a wide variety 
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of studies of galaxy formation and evolution (e.g., Croton 



et a l. 2006; B ower et ai:i[2006| |Somerville e t al. 2008; Dekel 
et al. 2009; Bouche et al.||2010|"]B'enson||2010 Hopkins et 



al. 2010| . They also complement basic analytic tools like 



the Press-Schechter and excursion set formalisms (|Press & 
Schechterl 19741 [Bond et al.|1991| ). 



Baryons, especially in gaseous form, are affected by a 
wealth of additional physical processes, including pressure 
forces, radiative heating and cooling, shocking, and star for- 
mation. It is therefore not surprising that they display a 
much richer phenomenology. Understanding the behavior of 
baryons is at least as fundamental as understanding the un- 
derlying dark matter, as they shape the galaxies that we 
observe. Further impetus for systematically quantifying how 
baryons assemble into dark matter halos is provided by re- 
cent work on the character of the accretion flows that fuel 
galaxies. Numerical simulations in fact support a picture in 
which galaxies acquire their gas through two distinct chan- 
nels, the "cold" and "hot" modes (Katz ct al. 2003 
et al.|2005[|Keres et al.|2009a||Birnboim fc D ckcl 2003] 



Keres 



Dekel 



fc Birnboim|2006| [Dekel et aL|2009| ). While this bimodality 
is closely related to the classical "rapid cooling" and "slow 
cooling" regimes (|Rees fc Ostriker|1977[|Silk|1977[ [White fc 
Rees[|1978 White & Frenk 19911, recent studies have shed 
new light on the properties of the flows and the transition 
between the two modes. 

According to the simulations, cold accretion (which we 
define as the mode in which the gas maintains a temperature 
T < 2.5 X 10^ K deep within the halo) is the main driver of 
galaxy formation in the sense that most of the baryons that 
make it into galaxies are accreted via this mode. At high 
redshift, the cold mode proceeds distinctly via filamentary 
streams. The existence of the cold mode and its importance 
in driving galaxy formation (in analogy with the rapid cool- 
ing regime) has been confirmed in analytic studies and in ID 
simulations ([Forcada-Miro fc White|1997[[Birnboim fc Dekel[ 



2003 Dekel fc Birnboim[[2006 l, in addition to cosmological 
simulations by different groups using both smooth parti- 
cle hydrodynamics (SPH) and grid-based codes (e.g. |Ocvirk[ 
et al.|2008; Brooks et al.|2009| [van de Voort et al.|2011[ ).lt 



is therefore a robust theoretical prediction. The hot mode, 
in which the accreting gas shock heats to roughly the virial 
temperature before cooling in a more spherically symmetric 
fashion onto galaxies, is predicted to compete with the cold 
mode only around the transition halo mass Mh ~ 3 x 10^^ 
M0 at 2 < 1 ( [Keres et al.|2009a l. 

As the cold streams penetrate the halos on a free fall 
time, they are very effective at transporting gas to their 
centers and may be connected to a variety of observed 
phenomena, including high-redshift star-forming galaxies 
(Elmegreen & Elmegreen 2006; Genzel et al. 2006; Dekel 
et al. 2009a,b), Lya blobs ( e.g.,|Fardal et al.|200T|[Dijkstra 
|fc Loeb|[2059l [Goerdt et al.||2010[ but see Faucher-Giguere 
et al. 2010), dense absorption systems in quasar and galaxy 



spectra (e.g., Dekel et al. 


2009 Faucher-Giguere fc Keres 


2011 Kimm et al. 2011 


Stewart et al. 20111, and high- 



quist|2009 \ . The morphological distinction between the cold 
and hot modes is also directly relevant to the efficiency of 
feedback processes in galaxy formation, as the physical re- 
quirements for preventing the infall of dense, cold streams 



or clumps are different from the requirements to prevent a 
spherical atmosphere from cooling. 

While previous work has investigated the detailed 
physics of gas accretion (e.g., [Keres et al.|2005[ [Keres et alTj 



2009a Ocvirk et al. 2008 1, with some broad considerations of 
the accretion rates into halos and galaxies, no previous study 
has systematically quantified the baryonic accretion rates a 
function of halo mass and redshift, for the different compo- 
nents, in a manner analogous to the growth and merger rates 
for the dark matter. As a consequence, theoretical analyses 
are often based on simplifying assumptions, such as assum- 
ing that the gaseous accretion or star formation rates follow 
the dark matter with a constant proportionality factor over 
some releva nt mass range (e.g., [Dekel et al. [2009 



Bouche 



let al."20 10l [Bauermeister et al.||2010[ [Krumholz fc BiirE" 
iert,|2010, jGattaneo et al.[[2010 1 . An important exception is 



the recent study of van de Voort et al. (2011 1, who quantr 



fied many of the relevant quantities using hydrodynamical 
simulations from the OWLS project ( [Schaye et al.|20To| in 
a general study of gas accretion onto galaxies and halos. 
Our focus here is however complementary in that we ex- 
tend the general studies of [Keres et al.[ ( [2005| ) and [Keres et[ 



al 



(2009a 



using a new set of higher resolution simulations 
and with different outflow parameterizations, to specifically 
quantify the accuracy of common approximations based on 
dark matter dynamics. Our study of halo accretion is also 



more detailed than that of van de Voort et al. ( 2011 1 in cer 



tain respects, as we explore the dependence on the reference 
radius, cover a wider redshift interval, and explicitly quan- 
tify the accuracy of several common scaling relations. Our 
simulations have a mass resolution better by a factor of four 



than the ones analyzed by Oppenheimer et al. (20101, who 



investigated stellar mass growth in hydrodynamical simula- 
tions but did not focus on the accretion rates. 

The richer physics of gas accretion, relative to the dark 
matter, introduces a complication: the accretion rates can be 
quite sensitive to feedback, especially in the form of galac- 
tic winds. Outflows must be considered in galaxy forma- 



et al.||2000||Pettini et al.||2001| [Shapley et al.||2003||Martin 


& Bouche[|2009 


Steidel et al.|[2010l, expected theoretically 


(e.g., Efstathiou 


2000 Murray et al. 2005 20111, and in- 



ferred from comparison of the measured star formation his- 
tory and stellar mass density with cosmological simulations 
(Springel fc Hernquist 2003; Keres et al. 2009b; Schaye et 
al. 2010; Oppenheimer et al. 2010). Unfortunately, the ex- 
act nature and properties of feedback remain major unsolved 
problems in galaxy formation, so that accretion rates mea- 
sured from hydrodynamical simulations risk being model de- 
pendent. In this paper, we demonstrate the dependence on 
the outflow parameters using a suite of otherwise identical 
simulations. We show that, even in the absence of outflows, 
the cold gas accretion rate is in general not a universal (or 
constant) fraction of the dark matter accretion rate. In the 
presence of outflows, the cold gas accretion rate itself can 
decouple from the star formation rate within the halo. In 
general, the star formation rate within low-mass halos is 
suppressed with respect to the cold gas accretion rate be- 
cause the gas can be ejected and unbound from the halo 
before forming stars. Additionally, the energy and momen- 
tum injected in the halo by outflows can in some regimes 
directly slow down gas accretion. Sufficiently massive halos 
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Table 1. Hydrodynamical Simulations 



Name 






e {/i-l kpc)': 


if 




?;w (km s 


gdm40n512 


40 


2 X 5123 


1.6 











gdm40n512_winds 


40 


2 X 512^ 


1.6 


1 


0.25 


342 


gdm40n5 1 2 _swinds 


40 


2 X 5123 


1.6 


2 


0.5 


342 


gdm40n512_fwinds 


40 


2 X 5123 


1.6 


2 


2 


684 



"Comoving box side length. 
''Total number of dark matter+gas particles in the box. 
'^Comoving Plummer equivalent gravitational softening length. 
'^Wind mass loading factor. 
■^Fraction of supernova kinetic energy injected into winds. 
•'Wind velocity (eq. (jijl). 

Note. — Dark matter, gas, and stellar particle masses are 3 X lO'^ Mq, 6 X lO'' h"^ Mq, and 3 X 10® Mq, respectively. 



can however retain the gas ejected from the galaxies they 
contain, and re-accretion can boost the star formation rates 



in those halos with respect to the no-feedback case (e.g., Op- 



penheimer et al. 2010|. In addition to the accretion rates. 



we also quantify the baryon mass fraction within halos, and 
show that it is also strongly suppressed by outflows in low- 
mass halos. Our simulation without outflows resolves the 
mass scale at which photoionization by the UV background 
suppresses the total baryon fraction at 2 < 5, corresponding 
to a circular velocity Vc ~ 50 km s~^. We provide simple fit- 
ting formulae for the cold gas accretion rate as a function of 
halo mass and redshift in the no-wind case, and compare it 
to recent star formation rate estimates of z ~ 2 galaxies. Our 
results should be useful to improve the accuracy of analytic 
(e.g.,|Boucheetal.|2010|, semi-analytic (e.g., |Cattaneo et al, 
" " "'~ ~ |Dutton et al. 



2006 



Khochfar & Silkl 



2009) |Lu et aL][20TT 



2010 [Benson fc Bower|20lT l, and empirical (e.g., Conroy & 
Wechslor 2009') models of galaxy formation. 

The plan of this paper is as follows. In i|2] we describe 
our numerical simulations and methodology. The accretion 
rates and baryon fractions measured from our suite of sim- 
ulations are presented in !|3j where we also discuss the im- 
plications of our results for dark matter-based scalings and 
provide relevant fitting formulae. We discuss our results and 
conclude in Q 



2 SIMULATIONS AND METHODOLOGY 
2.1 Code Details 

Our cosmological hydrodynamical simulations, run with the 
GADGET code (|Springel|2005[ ), are similar to the ones pre- 
sented in Keres et al. (2009a I and Faucher-Giguere et al. 
( |2010[ ) (the gdin40n512 simulation analyzed here is the same 
as in the latter article). Gravitational forces are calculated 
using a combination of the particle mesh algorithm (e.g., 
[Hockney fc Eastwood|[l988| for large separations and the 
hierarchical tree algorithm (e.g., Barnes fc Hut|1986 Hern 



quist 19871 at small distances. The gas dynamics is com- 



puted using a smoothed particle hydrodynamics (SPH) al- 
gorithm (e.g., |Lucy||1977||Gingold fc Monaghan||1977[ ) that 



conserves both energy and entropy ( [Springel fc Hernquist| 
20021. 



The modifications with respect to the public version 
of the code include the treatment of cooling, the effects of 
a uniform ultra-violet (UV) background, and a multiphase 



star formation algorithm as in Springel fc Hernquist ( 2003a I 



Star formation is implemented by the stochastic spawning 
of coUisionless star particles by the gas particles. In prac- 
tice, star formation in the multiphase model occurs above 
a density threshold of nn = 0.13 cm~3 and is calibrated to 
the observed |Kennicutt| ( |1998[ ) law. The thermal and ion- 
ization properties of the gas are calculated including all the 
relevant processes in a plasma with primordial abundances 



of hydrogen and helium follo wing |Katz et al. ( 1996 1. We use 
the UV background model of Faucher-Giguere et al. (20091, 



empirically calibrated to match IGM absorption measure- 
ments (Faucher-Giguere 2008a, b,c). 



2.2 Galactic Winds and Simulation Parameters 

To investigate the effects of feedback in the form of galac- 
tic winds on the baryonic accretion rates, we use a simple 
constant-velocity wind model, implemented as in ISpringel] 
|fc Hernquist I ( |2003a[ ). Exploring these effects is important, 
as we have outlined in the introduction, since outflows are 
known to be ubiquitous in the Universe (especially at high 
redshift), and to have important implications for the pre- 
dictions of numerical simulations. Our goal here is to show 
phenomenological examples of how galactic winds can affect 
the accretion rates of the different baryonic components. In 
particular, we do not imply that constant-velocity winds are 
an accurate representation of actual galactic winds, nor do 
we attempt to match the details of observations. Instead, we 
simply explore variations of the model parameters to deter- 
mine features that depend on (or are robust to) the wind 
properties. In future work, we will use more physical im- 
plementations of galactic outflows to study their dynamical 
interactions with accreting material in detail. 



As described in Springel fc Hernquist ( 2003a I , the wind 



model is parameterized by a constant mass loading factor rj 
and a constant fraction of supernova energy deposited into 
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Figure 1. Comparison of the median accretion rates through the virial shell for the different wind prescriptions. Solid curves: no winds. 
Dashed: constant-velocity winds with i)w = 342 km s~^ and mass loading rj = 1 (winds). Dotted: constant-velocity winds with i;w = 342 
km s~^ and mass loading rj = 2 (swinds). Dash-dotted: const ant- velocity winds with Uw = 684 km s~^ and mass loading rj = 2 (fwinds). 
The thick grey lines show the fit of FMBIO to the dark matter growth rate, scaled by (Hm — f!j3)/f2m since these authors studied a dark 
matter-only simulation. 



the wind x- Specifically, the disk mass-loss rate in the wind 
is defined as — rjMf,, where = SFR is the star 
formation rate, and the velocity of the wind as it leaves the 
disc satisfies the equation 
1 



I.e., 



2X£SN 

n 



(1) 



(2) 



In the equations above, esN ~ 4.6 x 10*** erg Mq^ is the aver- 
age energy output by supernovae per unit stellar mass, corre- 
sponding to the canonical value of 10^^ erg for a Salpeter ini- 
tial mass function (IMF) between 0.1 and 40 M©, assuming 
that stars with > 8 Mq die as supernovae. Gas particles are 
stochastically incorporated into the wind according to their 
star formation rate and the wind particles are briefly de- 
coupled from hydrodynamic interactions. This approach has 
been shown to yield numerically converged results even at 
the relatively coarse resolution of cosmological simulations 
(Springel & Hernquist"2003a|. In order to produce roughly 



bipolar outflows normal to galactic discs, the wind particles 
are given initial velocities parallel to ±V(^ x v, where (j) is 
the local gravitational potential and v is the velocity of the 
star-forming particle sourcing the wind. 

Our simulations implement a basic form of instanta- 
neous gas recycling from Type II supernovae, in which a 



fraction /3sn = 0.1 of the star formation rate is immediately 
returned as interstellar gas. Simulated galaxies can therefore 
maintain steady-state star formation rates exceeding the ex- 
ternal gas supply rate by a fraction Psn- The additional star 
formation however does not contribute to the net rate of 
stellar mass growth and the parameter Psn is sensitive to 
uncertain IMF and stellar evolution assumptions. In order 
to subtract out the model dependence on /3sn, we henceforth 
make the identification SFR S'F7?'°*(1-/3sn) unless oth- 
erwise specified. I.e., the SFRs we report do not include the 
boost from gas return. This convention has the benefit that 
SFRs exceeding the cosmological supply rate are directly 
indicative of additional supply channels, e.g. via wind fall- 
back (either on the original galaxy, or on a neighbor) or hot 
gas cooling. To be exact, our SFRs should however be mul- 
tiplied by (1 — /3sn)~^ ~ 1.1 when comparing with measured 
star formation rates. 

The simulations analyzed in this work are detailed in 
Table [T] The cosmological parameters for all these simula- 
tions are the be st-fit parameters for the WMAP5-I-BAO-I-SN 
data ( [Komatsu et al._ 2009J : (Jim, ^h, ^a, h, as, ris) = 
(0.28, 0.046, 0.72, 0.7, 0.82, 0.96). All simulations have a 
box side length of 40 comoving Mpc and use 512^ par- 
ticles to represent each of the dark matter and gas compo- 
nents. They assume alEi senstein fc Hu| ( |1999| l transfer func- 
tion and the A''— body integration begins at redshift z = 99. 
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The dark matter, gas, and stellar particle masses are 3 x 10 
h-^ Mq, 6 X 10*^ h-^ Mq, and 3 x lO'' Mq, respectively. 

The basic simulation gdm40n512 does not include feed- 
back, except for a sub-resolution multiphase model to pres- 



surize the interstellar medium (ISM) ( Springel & Hernquist 
2003a I. Three otherwise identical simulations incorporate 



galactic winds as described above, with parameters given in 
the table. Briefly, the winds in the winds and swinds runs 
have a launch velocity of «w ~ 342 km s~^, and mass loading 
factors rj — 1 and rj — 2, respectively. The swinds model is 
similar to the slow winds (sw) of Oppenheimer et al. (20101, 
which they find produce agreement with the observed stellar 
mass fraction at z = comparable to their best momentum- 
driven winds (though not as good at the low-mass end).[^ 
The third wind simulation (fwinds) uses an extreme set of 
parameters, rj = 2 and x = 2, which corresponds to driving 
the winds with twice as much kinetic energy as supernovae 
fiducially provide, and ejecting twice as much mass in the 
winds as the star formation rate. The wind velocity in this 
case is 684 km s~^. Although other feedback mechanisms, for 
example radiation pressure from massive stars (e.g., Murray 




et al. 20051, can accelerate winds with velocities or mass 



loadings greater than possible for energy injection by super- 
novae alone, we have verified that these parameters yield a 
baryon mass function at z = heavily suppressed relative to 
observations (e.g., Bell et aL]|2003 1 for 2 decades in galaxy 
mass around L* (similar to the cw model of Oppenheimer 
et al. 2010). 

The simulations presented in this work have a mass res- 
olution 11 X finer th an the main simulation analyzed by 
3res et al. (2009a I, in a comparable volume (2 x 288'^ 
particles with 50 comoving Mpc box side length vs. 
2 x 512^ particles with a 40 comoving Mpc box side 
length), and use the same code. Our finer mass resolution 
is also slightly better than that of the smaller 11 co- 



10.5 11.0 11.5 12.0 12.5 13.0 13.5 



Figure 2. Comparison of the cold gas accretion rates through 
the virial shell for the different wind prescriptions. Only the gas 
particles with inward radial motion are included in order to ex- 
clude the negative contributions from outflowing material. From 
the bottom up, the results are shown at z = 0, 2,4, and 6. Solid 
curves: no winds. Dashed: constant- velocity winds with v-w = 342 
km s~^ and mass loading rj = 1 (winds). Dotted: constant- velocity 
winds with = 342 km and mass loading rj = 2 (swinds). 
Dash-dotted: constant- velocity winds with Dw = 684 km s~^ and 
mass loading rj = 2 (fwinds). 



our results. As our simulations resolve the UV background 



moving Mpc box (with 2 x 128 particles) used by|Keres et] suppression scale at redshifts z < 5 (p^, there is only a 



al. (2009a I to test the numerical convergence of their re- 
sults, so that we can be confident that our findings, broadly, 
are robust. It is worth noting that many of the quantities 
that we analyze, being evaluated at the virial radius and 
therefore involving lower densities, should actually have sig- 
nificantly less stringent convergence requirements than the 



accretion rates onto galaxies studied by Keres et al. ( 2005 1 



and |Keres et al.| ( |2009a[ ). At the detailed quantitative level, 
it is however possible that even the accretion rates through 
the virial shell could be subject to minor resolution effects. 
Two effects stand out in particular. First, a higher resolu- 
tion allows the formation of lower-mass objects, which can 
remove gas by locking baryons into stars. Because structure 
formation proceeds hierarchically, the effect can propagate 



small regime of halo mass and redshift for which increasing 
the resolution might yield more star formation in low-mass 
halos. Furthermore, at z ^ 5 and for the minimum halo 
mass Mh ~ 10^'''^ M© at which we report accretion rates, 
the accretion rate of stars in our no- wind simulation is only 
< 3% of the cold gas accretion rate. The convergence re- 
quirements for the simulations with outflows are less severe 
in this respect because outflows are effective at removing 
baryons from low-mass halos. A second potential effect con- 
cerns high-mass halos, in which gas from a hot atmosphere 
can cool onto galaxies and also fuel star formation. This may 
in fact be somewhat sensitive to resolution, but in this case 
the situation is less clear and different numerical methods 



to higher-mass halos. Keres et al. (2009a I saw this effect 



at the 30% level on the smooth gas accretion rates. This 
effect should be significantly smaller in our simulations for 
several reasons. Since we do not explicitly distinguish be- 
tween smooth and clumpy accretion, it is only the fraction 
of baryons locked into stars in low-mass halos that can affect 



have not yet converged to a common answer (e.g., Keres et 
[al]|2009a[ ). This issue will be explored by comparing GAD- 
GET results with the moving-mesh code AREPO in a sub- 
sequent paper (Keres et al., in prep.). 



2.3 Halo Identification 



^ The simulations of [Oppenheimer et al.| | |2010| differ slightly 
by their inclusion of metal cooling and asymptotic giant branch 
(AGB) gas recycling, but we have verified that our swinds simu- 
lation produces a consistent mass function at z = 0. 



A friends-of-friends (FoF) algorithm (e.g., Davis et al.|1985 1 
with linking length set to b — 0.2 in units of the mean 
interparticle separation is used to identify the dark mat- 
ter halos in the simulations. The total mass of the par- 
ticles within each FoF group, Mpof, corresponds approx- 
imately to the mass in a sphere of mean interior density 
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Figure 3. Net cold gas accretion rate as a function of redshift 
and halo mass, for the simulation without outflows, through shells 
of radii l-Rvin O.S-Rvin ^.nd Q.2Ryi^. The panels show how deep 
the cold gas penetrates into halos in different regimes. The ver- 
tical axes cover exactly 2 orders of magnitude in each panel, so 
that relative vertical separations in the different panels are fair 
representations of ratios. 



180 times the background matter density, Mi^ob- |White| 



( 2002 1 in fact showed that Migot is a more accurate proxy 



for the FoF mass than the common virial mass definition 
of 200 times the critical density. We therefore define the 
virial radius of a halo as the radius of a sphere contain- 
ing Mi80b, Rvir = Risob ~ [MFoF/2407rp„(^)]^''^ where 
pu(z) — pcritf2m(l + z)'^ aud Pcrit is the critical density 
at z = 0. In what follows, we use Mh as a shorthand for 



MpoP ~ Migob- The center of each halo is defined as the 
point deepest in the gravitational potential. Since we run 
the FoF algorithm on the dark matter particles only, we 
multiply the returned masses by nm/(nm — fib) to account 
for the baryons. The use of the FoF mass is consistent with 



the previous work of McBride et al. ( 2009 1 and Fakhouri 



et al. (20101 on the dark matter mass assembly, and also 



allows for the direct use of existing analytic approximations 
to the halo mass function (such as that of Seth & Tormen 
2002, which is calibrated to FoF masses) in applications. 



2.4 Definition of Accretion Rates 

Our approach to calculating the accretion rates differs from 
the technique often used in Lagrangian simulations (either 
dark matter-only or SPH) of identifying particles belonging 
to halos or galaxies at different redshifts and dividing the 



mass difference by the time interval (e.g., McBride et al. 
2009; Fa khouri et al.||2010| [Murah et al.||2002| [Mailer et al. 
2006[|Keres et a l.'2009a', |van de Voort et al.|2011| ). Since the 
thermal state of the gas, in particular whether it is classified 
as "hot" or "cold" , can be a strong function of radius from 
the halo center, it is advantageous to measure the accretion 
rates through shells of prescribed radii. This allows us to 
investigate how the mix of different baryonic components 
evolves during the infall. Furthermore, this approach allows 
us to select particles based on the direction of their motion, 
which as we will see is useful for separating accreting mate- 
rial from outflows. Mass fluxes also have the significant ad- 
vantage that they can be identically computed in grid-based 
simulations (e.g., Ocvirk et al.|200"8 l; results from this work 
can therefore be directly compared with those obtained with 
grid codes, including the new moving-mesh code AREPO 
(Springel 2010 ). Numerical noise from the finite number of 
particles overlapping with a relatively thin shell is mitigated 
by evaluating the accretion rates only for halos containing 
more than 500 dark matter particles and, on average, an 
equal number of gas particles. Furthermore, all of our re- 
sults regarding low-mass halos are aggregate statistics over 
large samples and so are robust. In this work, we do not ex- 
plicitly distinguish between mergers and smooth accretion, 
but we do consider gas that is part of the ISM separately 
from the rest (see below), which we subsequently refer to as 
"diffuse", as explained below. 

Since the natural scale at which to evaluate the ac- 
cretion rates is the virial radius, we will usually refer to 
the "virial shell", although in general we can set its radius 
Rs{M, z) to any factor of the virial radius. The thickness of 
the shell is labeled by A_Rs(M, z), and it is centered on the 
point deepest in the halo potential. The net accretion rate 
of a given component i is then defined as 



VA 



fh - rp 

|rh - rp 



(3) 



where the sum is over the particles of the component of inter- 
est that intersect the virial shell, i.e. Rs — ^ |rp; — rh| ^ 



R, + 



Here, Mp; is the particle mass, rp; is its position 



vector, Vpj is its velocity vector, and rh is the position of 
the halo center. The velocity vector va, to be discussed in 
detail in the next section, accounts for the motion of the 
area element on virial shell itself. The velocity vector Vpj in- 
cludes both the peculiar velocity contribution and a radial 
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component H{z)Rs (assumed constant within a given radial 
shell) accounting for the Hubble expansion. Unless otherwise 
noted, all quantities are in proper units. Note that the con- 
tribution of any particle to the net accretion rate is positive 
if it is moving inward with respect to the virial shell, and 
negative otherwise. We will also consider the total mass en- 
closed within the virial shell of different components, defined 
as 



(4) 



where the sum is over all particles with | rpj — j ^ -Rs . The 
star formation rates within halos are evaluated similarly. 

We keep track of five species separately: dark matter 
(DM), cold gas (nn < 0.13 cm"^, T sC 2.5 x lO'^ K), hot 
gas (riH < 0.13 cm"^, T > 2.5 x 10^ K), stars (★), and ISM 
(riH 0.13 cm~^). The temperature threshold between the 
cold and hot phases is motivated by the bimodality found 
in cosmological simulations for the maximum temperature 



attained by gas accreting from the IGM (e.g., Keres et al 



20051, and approximately distinguishes between gas that is 
shock- heated as it falls in, and gas that is not significantly 
shocked. Note, however, that the cold and hot gas accretion 
rates through shells considered in the present work differ 
from the "cold mode" and "hot mode" defined by | Keres | 
et al. (20051 and Keres et al. (2009a I. In those papers, the 



cold and hot modes are defined in terms of the maximum 
temperature attained before accretion onto galaxies, rather 
than the instantaneous temperature. The ISM material is 



defined in the context of the jSpringel & Hernquist ( 2003a I 



sub-resolution multiphase star formation model used in our 
simulations. As outlined in |2.1[ the gas particles with den- 
sity nil J? 0.13 cm~^ are assumed to develop a multiphase 
structure. These multiphase particles carry an effective tem- 
perature that is a mass- weighted average of cold and hot ISM 
phases. As a result of the averaging procedure, these effec- 
tive temperatures can exceed the T = 2.5 x 10^ K threshold 
between the cold and hot phases, even though most of the 
gas mass is usually cold. We therefore treat the ISM material 
separately and never include it in the cold and hot compo- 
nents, which therefore always have have densities nu < 0.13 
cm~'^. Like cold accretion from the IGM (but unlike hot gas 
in massive halos) , gas that is accreted by a halo in ISM form 
(e.g., via mergers) can rapidly form stars. It will therefore 
be useful to consider the sum of the cold and ISM accretion 
rates in some circumstances. The net accretion rates are fur- 
ther divided into "inward" and "outward" parts, such that 
Mi = Mi^ - M°"*. The inward rates are defined by sum- 



> 0, as 



ming only over particles with (vpj — va 
in equation (j3|, and vice versa. 

Because our Eulerian definitions of the cold and hot gas 
accretion rates differ from the Lagrangian definitions of the 
cold and hot modes of Keres et al.| ([2005 J and [Keres et al 



(2009a I, some systematic differences are anticipated. Since 



we have adopted the same temperature threshold between 
the cold and hot components as in these previous papers, 
but consider the instantaneous instead of maximum tem- 
perature, we expect our cold gas accretion rates onto halos 
to be systematically higher. A similar consideration applies 
when comparing with the Lagrangian definition of [van de] 
[Voort et al.| ( |2011[ ), but some deviations from this expec- 
tation could arise due to the higher temperature threshold 



T = 10^'* K adopted by those authors. Where our results 
overlap, they are broadly consistent with those previous 
studies and with the adaptive mesh refinement simulation 



of Ocvirk et al. (20081, as we will detail where applicable. 



We note, however, that many of our results are new and not 
directly comparable with earlier SPH analyses. For exam- 



ple, Keres et al. (2009a I focused on accretion onto galaxies 



rather than halos, and van de Voort et al. (20111 did not 
resolve the accretion as a function of radius or based on 
whether the gas elements were inflowing or outflowing. 

It is also noteworthy that the simulations of |van de| 
Voort et al. (20111 implemented some physics differently 



from our simulations. In particular, their reference model 
included metal line cooling, which they found to have an 
important effect on accretion onto galaxies, but a much 
smaller impact on halo accretion. An importance distinc- 
tion between the cold vs. hot accretion fractions reported 



herein and in Keres et al. ( 2009a I is that in the latter work 



only gas that ended up into galaxies was considered. Since 
the gas accreted by galaxies is strongly biased toward the 
cold mode, especially in massive halos in which the cool- 
ing times for gas shock-heated to the virial temperature are 
long, the fraction of cold vs. hot mode accretion onto galax- 
ies can be much larger than the relative cold vs. mode hot 



rates through the virial radius. As we will show (i 3.2.2 1, 



the cold gas accretion rates at the virial radius is closely re- 
lated to the star formation rates within halos, and therefore 
share many of the key features of the Lagrangian cold mode 



defined by Keres et al. ( 2005 1 and Keres et al. ( 2009a I 



Since our accretion rates are evaluated directly using 
the velocity vectors of the particles, they are instantaneous. 
On the other hand, accretion rates evaluated by dividing 
the mass by the time interval between discrete snapshots are 
implicitly time-averaged. While time averaging preserves the 
mean accretion rates, it in general modifies the shape of the 
distribution, reducing its scatter, and can therefore affect 
statistics such as its median. This effect will be encountered 
in §3.1| where we compare the median dark matter accretion 



rates inferred from the Millennium simulations to ours. 

Plotting the halo accretion rates and baryon contents 
as a function of halo mass requires some care, since the slope 
of the halo mass function implies that the data points in any 
given mass bin are in general concentrated in its lower end. 
To mitigate this potential source of bias, we use the median 
Mh in each bin as the horizontal coordinate, rather than its 
geometric center. 



2.5 Definition of the Virial Shell 

The virial shell in general grows with the Hubble expansion 
and as the halo acquires mass, and the shell also moves with 
the halo itself. An area element on the shell therefore has 
motion of its own, which must be accounted for in calculat- 
ing fluxes. The velocity vector of a shell area element can be 
written as 



(5) 



where the Vh term accounts for the center-of-mass motion of 
the halo and the Vcxp term accounts for the expansion of the 
virial shell. Let ta be the position vector of the area element 
of interest on the virial shell (assumed to be at the center 
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Figure 4. Comparison of the ratio of the cold gas+ISM to dark matter accretion rates through the virial shell for the different wind 
prescriptions (blue). The red curves show the same ratio, but including all gas. Solid curves: no winds. Dashed: constant- velocity winds 
with = 342 km and mass loading rj = 1 (winds). Dotted: constant- velocity winds with v^, = 342 km and mass loading rj = 2 
(swinds). Dash-dotted: constant- velocity winds with v^j = 684 km and mass loading r) = 2 (f winds). The thick grey lines show the 
universal ratio Qb/f^DM- 



of its thickness). Then, since the shell expansion velocity is 
proportional to the rate of change of the shell radius, 



FA - rii 
|rA - Th 



(6) 



Although the virial shell can be defined arbitrarily, for def- 
initions tied to the virial radius of halos its total rate of 
change is a sum of the rate of change arising from the Hub- 
ble expansion and a term arising from the growth of the 
mass in the halo: 



■ _ dFU dAh dRsdz 
' ^ 9Mh dt Ihdt' 



(7) 



This derivation is completely general as long as Rs = 
7?s(Afh, z), and any redshift dependence of the density con- 
trast used to define Rs is folded in the second term on the 
right-hand side. For R^ — Rvii = Risob, straightforward al- 
gebra yields 



Risob = -Ri 



Hiz) 



Misob 



3A/l: 



(8) 



We can use the fitting formula of Fakhouri et al. (20101 
(FMBIO) for the average mass accretion rate of halos mea- 
sured from a dark matter-only simulation to estimate the rel- 
ative magnitude of the two terms in this expression. FMBIO 



found that the median mass accretion rate follows 

Mvl^j (1 + 1-65.) 



(Afi80f,)mod ~ 25.3 Mq yr-^ 



(9) 



where we use the « sign to indicate the fact that FMBIO 
actually measured the growth of the FoF mass, and not di- 
rectly the growth of Misot., although the two closely follow 
each other, as we have emphasized above. Using this equa- 
tion gives 

0.1 

0.12 1 ,::;t: \ (i + i.es.), (lo) 



mcd 



3Mi80bH{z) 



JVllSOb 

10^2 M<7 



implying that the Hubble expansion term dominates at low 
z and for Migoft ^ 10^^ M©, but not necessarily at higher 
mass and redshift. 

To compare with previous results on halo mass growth, 
such as those of FMBIO, the most relevant definition is a La- 
grangian shell that follows the halo and grows with it. We 
therefore include both the Hubble expansion and the mass 
growth term in modeling the expanding virial shells. The 
mass growth term is calculated on a halo-by-halo basis, us- 
ing the mass growth rate measured for each halo. Unless oth- 
erwise specified, we adopt -Rs — Rvir and A.Rs = 0.2A_Rvir 
throughout, which we have found to be a good compromise 
between numerical noise and spatial resolution. 



© 0000 RAS, MNRAS 000, 000-000 



z=0 



z=1 



Baryonic Halo Assembly 9 
z=2 




z=3 



z=4 



z=5 



- * 




i 






t' 







10.5 11.0 11.5 12.0 12.5 13.0 10.5 11.0 11.5 12.0 12.5 13.0 10.5 11.0 11.5 12.0 12.5 13.0 

logio M, (M g) logio M, (M g) logio M, (M a) 



Figure 5. Comparison of tlie ratio of the star formation to cold gas+ISM accretion rates through the virial shell for the different wind 
prescriptions. Solid curves: no winds. Dashed: constant-velocity winds with Hw = 342 km s~^ and mass loading r] = 1 (winds). Dotted: 
const ant- velocity winds with v^i = 342 km s~^ and mass loading rj = 2 (swinds). Dash-dotted: constant-velocity winds with Hw = 684 
km s"-"^ and mass loading r) = 2 (f winds). The thick grey lines show the reference value of unity. 



3 RESULTS 



3.1 Accretion Rates 



We begin by quantifying the accretion rates of the dark mat- 
ter, cold gas (riH < 0.13 cm"^, T < 2.5 x 10^ K), hot gas 
(nn < 0.13 cm~^ T > 2.5 x 10^ K), and ISM (nn ^ 0.13 
cm~'^) as a function of halo mass and redshift. For compari- 
son, we also present results for the star formation rate within 
the halos in the simulations. Figure [T] shows the median net 
accretion rates through the virial shell, for the simulations 
with different outflow parameterizations, at z = — 5. 

As expected, the dark matter accretion rates are insen- 
sitive to the wind feedback, since they are not subject to 
hydrodynamical forces. The median net dark matter accre- 
tion rates from our hydrodynamic simulations are in excel- 
lent agreement with the fitting formula of FMBIO (eq. (|9|) 
based on the Millennium simulations, at all redshifts 2^1. 
Differences between the median dark matter accretion rates 
of a factor « 1.5 — 2 are apparent at z = 0. However, our 
mean dark matter accretion rates are in excellent agreement 
with the corresponding FMBIO fit at all redshifts, indicat- 
ing (as we have verified) that it is the shape of the dark 
matter accretion rate distribution that is slightly different 
in our simulations. The dominant effect likely lies in the 
definition of the accretion rates. In this work, the accretion 
rates are instantaneous and directly evaluated from particle 
velocities in single snapshots. FMBIO, on the other hand, 
evaluated the accretion rates by dividing the mass differ- 



ence in subsequent snapshots by the time interval. At « = 0, 
the Millennium snapshots are spaced by ~ 260 Myr and 
so FMBlO's accretion rates are averages over a comparable 
time scale. The effects of averaging are most pronounced at 
z ~ 0, where accretion events are rarest, explaining why 
they are not as apparent at z ^ 1. Additional modest differ- 
ences are expected due to the different (first-year WMAP) 
cosmology assumed in the Millennium simulations. 

The baryonic rates, on the other hand, depend substan- 
tially on the wind prescription. In particular, the star forma- 
tion and net cold gas accretion rates in low-mass halos are 
strongly suppressed as winds remove gas from the galax- 
ies. In sufficiently massive halos, however, Figure [l] shows 
that the accretion rates of the cold and hot gas are largely 
insensitive to the wind models implemented in the simu- 
lations. This independence arises in part because as halo 
escape velocity becomes sufficiently large, the outfiowing 
material cannot reach the virial radius. As demonstrated 



Oppenheimer & Dave ( 2008 1 and Oppenheimer et al 



(20101, hydrodynamic interactions with the gas surrounding 



the galaxy also help confining outflowing material in high- 
mass halos. These authors further showed that, as a result, 
massive halos preferentially confine their outfiows even in 
models in which outfiow velocity scales with the escape ve- 
locity of the halo, as expected in momentum-driven models 
(e.g., [Murray et al.|2005 |. 

It is important to note that although the gas accre- 
tion rates through the virial radius appear insensitive to 
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Figure 6. Comparison of tiie ratio of the star formation to dark matter accretion rates through the virial shell for the different wind 
prescriptions. Solid curves: no winds. Dashed: constant-velocity winds with Hw = 342 km s~^ and mass loading r] = 1 (winds). Dotted: 
const ant- velocity winds with v^i = 342 km and mass loading rj = 2 (swinds). Dash-dotted: constant-velocity winds with Hw = 684 
km s"-"^ and mass loading r) = 2 (f winds). The thick grey lines show the universal ratio f2b/f^DM- 



the details of the outflow prescription in massive halos, star 
formation feedback has been typically found to be insuf- 
ficient in suppressing the galaxy baryon mass function to 
the observed level in these galaxies (e.g., Oppenheimer et 
al.|2010 1. Further feedback by active galactic nuclei (AGN), 



not included here, is probably responsible for suppressing 



the galaxy mass function in this mass regime (e.g., Croton 



et al.|2006||Somerville et al.||2008[ [Hopkins et al.|2008 i" 

Comparing the net and pure infall baryonic accretion 
rates is useful for assessing whether the net influx is affected 
by galactic winds solely because of the negative contribution 
from outflowing material, or whether infalling material can 
also be effectively kept out of halos via interactions with 
the outflows. If galactic winds simply eject material from 
the halos but do not signiflcantly disrupt infalling baryonic 
structures, then the pure infall rates should be insensitive to 
the wind prescription, even if the net accretion rates depend 
strongly on it. This comparison is particularly instructive for 
the cold gas component, as it is predicted to provide most 



of the fuel fo r star formation (e.g., Keres et al.|2005 Keres 



et al. |2009a| . As Figure [i] shows, the pure infall Af^^u rates 



are indeed less sensitive to the wind prescription than the 
net rates (c.f.. Fig. [T]), and relatively weakly dependent on 
it for the winds and swinds simulations, with Ww = 342 km 
s~^. For our more extreme swinds simulation with u„ — 684 
km s^^ and rj = 2, however, even the pure infall cold gas 
accretion rates are substantially affected at most redshifts, 
indicating that the outflows can in fact affect the infalling 



cold filaments. This is consistent with the inference of |Op-| 
ipenheimer et al. ( 2010 1 based on the suppression of the star 



formation rates by outflows (see also f 3.2.2 1, but is for the 
flrst time directly shown here. 

In Figure[3] we show how the net cold gas accretion rate 
varies with redshift and halo mass, for the no-wind simula- 
tions, but through shells of radii — i?vir, 0.5i?vir, and 
0.27?vir {ARs = 0.27lvir, 0.2i?vir, and O.liJvir). This illus- 
trates how deep the cold gas penetrates into halos in differ- 
ent regimes. While at z 2 the cold gas accretion rate is 
signiflcant down to 0.2i?vir even in halos of mass Mh > 10^^'^ 
Mq that are dominated by hot gas (^3.31, the cold gas ac- 



cretion rate at small radii is markedly suppressed in those 
halos at 2 5C 1. This conflrms the persistence of cold streams 
into massive halos at high redshift ( Keres et al.|2005 Keres 
et al.|2009a||Ocvirk et al.|2008| ), arising from the short cool- 



ing times in the dense gas brought in along the dark matter 



fllaments into rare halos ( Dekel & Birnboim 2006 1 . This phe- 



nomenon operates increasingly less efficiently at 2: < 2 as the 
densities become lower, the cooling times correspondingly 
longer, and halos of a given mass become more abundant 
and tend to reside inside large dark matter filaments, rather 
than at their intersection ( Katz et al.|[20 03). In addition to 
shocking, the decreasing cold gas accretion rates with de- 
creasing radius are caused by accretion by satellite galaxies 



on the way down to the center of the main halo (e.g., Keres 
et al.|2009a1|Simha et al.|2009[ ). 



In some regimes, particularly in high-mass halos, the 
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hot gas accretion rates are found to be insensitive to the 
outflow prescription (see also van de Voortet al.|2011[ ). This 
indicates that outflows have a small impact on gas shocking 
at the virial radius. However, this does not hold universally 
and Figure [l] shows significant outfiow dependences at lower 
masses, where winds can reach the virial radius with signif- 
icant velocities. 



3.2 Validity of Simple Scaling Relations 

One motivation for quantifying the baryonic accretion rates 
onto halos is to provide more realistic ingredients for galaxy 
formation models, i.e. to extend the standard framework 
based on dark matter-only results that has been the back- 
bone of most theoretical models so far. In this section, we 
assess the validity of several scaling relations that are of- 
ten used to model gas accretion and star formation in halos 
based on dark matter-only models. Such scaling relations 
have been qualitatively motivated by the results of hydro- 
dynamical simulations on the cold mode. Although limita- 
tions to such assumptions have already been clearly outlined 
(for example, hot mode accretion is more important in some 
regimes; e.g., Keres et al.|2005" Keres et al.|2009a l, the level 
of their accuracy has not been systematically quantified over 
broad intervals of halo mass and redshift, and in the pres- 
ence of outflows. We first review the motivations for such 
scaling relations, before exploring them in detail. 

Because the cold mode is characterized by infall that 
follows the dark matter filaments closely at large radii, it 
is tempting to assume cold gas enters halos at a rate fol- 
lowing the dark matter with a universal baryonic mass fac- 
tor, i.e. Mcoid + MisM ~ (nb/f^DM)AfDM. Furthermore, 
the cold mode is predicted to provide the main fuel for 
star formation. One might therefore assume that the star 
formation rates within halos are approximately equal to 
the rates at which they are supplied with cold gas, i.e., 
SFR « Mcoid+AfisM, as is the case globally in some regimes 



other galaxy. The (Mcoid + AfisM)/MDM ratio is strongly re- 
duced at halo masses Mh > 10^^'^ M0, as halos above this 



in the absence of feedback (e.g., Keres et al.|2005 Keres et 
aL]|2009a |. Together, these scalings lead to the common as- 
sumption that the star formation rates within halos is ap- 
proximated by a universal baryonic factor of the rate at 
which they acquire dark matter: SFR « (r2b/f^DM)AfDM 
(e.g., IGenel et al.||2008l |Dekel et al.||2009| |Bouche et al. 



2010), with an optional efficiency factor. It is however clear 
that these simple scalings cannot apply universally, as the 
cold mode is suppressed in halos above the transition halo 



mass Mt ~ 3 X lO" Mq (e.g., 


Birnboim & Dekel 2003 


Keres 


et al.|2005 Keres et al.|2009a 


Ocvirk et al.|2008 1. As shown 



above, feedback can also strongly affect both the gaseous ac- 
cretion and star formation rates. It is therefore important to 
quantify the applicability of these simplified assumptions. 



3.3.1 When Does Mcoid + Mism « {^b/^DM)MDM ? 

Figure |4] compares the cold gas+ISM and dark matter ac- 
cretion rates in our different simulations with and without 
galactic winds. While the gas accretion rates shown are dom- 
inated by the diffuse cold component, we have included the 
ISM contribution since most of it is also in cold form, and 
it can in principle be rapidly converted into stars as it falls 
into halos, either spontaneously or after merging with an- 



mass maintain stable virial shocks (e.g., p3irnboim fc Dekel 
2003}|Keres et al.|2005| [Keres et al.|2009a[ ) and" an mcreasmg 
fraction of the infalling gas is therefore hot, as indicated by 
the red curves which include all the gas. In high-mass halos 
where the total gas to dark matter accretion rate ratio is 
significantly below universal, the accretion of stars make up 
most of the difference. 

At lower halo masses, (Mcoid + Mi3m)/Mom is further 
suppressed in simulations with outflows. As discussed in §3.1 
(c.f., Figs[l]and[2]|, this is predominantly due to cold outflow- 
ing material contributing negatively to the virial shell flux, 
although hydrodynamical interactions between the inflow- 
ing and outflowing gas do suppress the infall rate of cold gas 
non-negligibly. At late times, Figure[4]shows that the median 
(A/cold + Mism) / Mom is systematically below the universal 
value Slb/f^DM even for low-mass halos that are not subject 
to strong virial shocking and in the absence of winds. In- 
cluding the hot gas does not significantly change this ratio, 
and we have verified that further adding stars (i.e., including 
all the baryons) does not change the picture at the low-mass 
end, so that more than half of the halos in this regime accrete 
baryons at sub-universal rates. In the absence of outfiows, 
the main mechanism keeping baryons out of low-mass ha- 



los is photoheating by the UV background (e.g., |Efstathiou 



T992{ [Thoul fc Weinberg|[T996l [Gnedin|[2000l) . This effect is 



quantified in §3.3| for our simulations, and is in fact sufficient 
to suppress the net accretion rate of baryons in halos with 
< 10^"-^ Mq at 2 sC 1. 

Interestingly, the z — 3 ~ 5 panels of Figure [4] show 
that halos with Mh < 10^^ Mq can accrete cold gas at 
"super-universal" rates with respect to the dark matter, i.e. 
(Mcoid + Mism) > r^b/i^DM, for the simulation without out- 
fiows. The UV background is also likely responsible for this 
effect: the gas kept out of low-mass halos remains available 
for future accretion and can contribute to super-universal 
gas accretion rates in other halos. Net accretion of baryons 
rates > (f2b/^^DM)MDM in some regimes is in fact required 
to explain the super-universal median baryon fractions over 
important intervals in halo mass and at all redshifts studied 
here, in the absence of outflows (j ]3.3[ ). This effect, albeit 
small and only apparent in simulations without strong feed- 
back, is also present in a number of other cosmological sim- 
ulations using different codes ( [Springel fc Hernquist|[2002| 
Keres|[2"007 \ and is therefore believed to be real. The tran- 



sition versus redshift from super- to sub-universal median 
baryonic accretion rates in the low-mass bin Mb ~ 10^"'^ 
M0 can be explained by the redshift dependence of the UV 
background suppression scale discussed in §3.3| In models 
with outflows, the ejection of baryons from low-mass halos 
could also in principle permit subsequent super-universal ac- 
cretion onto other halos, but this is not apparent in Figure 
[4] indicating that it is subdominant to the first-order effect 
of winds ejecting gas from the halos. 



3.2.2 When Does SFR « (Mcoid + Msm) ? 

Without strong feedback, star formation in galaxies is ex- 
pected to follow the cosmological supply of gas. Observa- 
tionally, galaxies otherwise exhaust their fuel on time scales 
~ 1 — 2 Gyr locally, with some evidence of even shorter 
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Figure 7. Comparison of tlie median baryon mass fractions within halos, broken down by components, for the different wind prescriptions. 
We comment on the relative contributions of stellar and ISM material in §3.3| Solid curves: no winds. Dashed: constant- velocity winds 
with i)w = 342 km and mass loading rj = 1 (winds). Dotted: constant-velocity winds with Dw = 342 km and mass loading rj = 2 
(swinds). Dash-dotted: constant- velocity winds with tiw = 684 km and mass loading r) = 2 (f winds). The thick grey lines show the 
universal ratio f7^,/nni. 



gas-depletion time scales ~ 0.5 Gyr at 2: ~ 2 (e.g., Genzel 
|et al.|2010[ ). Simulations adopting star formation laws moti- 
vated by observations in fact confirm that the star formation 
rates are limited by the accretion of intergalactic gas, and 
have further predicted that cold mode accretion is the main 
channel by which galaxies get their gas (e.g 



2002 Keres et al. 2005 Keres et al. 2009a I. This is illus- 



Murali et al. 



trated in Figure [5j in which we plot the median ratio of the 
star formation rates within halos to the rates at which cold 
gas-l-ISM enter their virial radii, SFR/{Mcoid + A/ism). At 
redshifts z ~ 1 — 3, the median ratio is close to unity at 
all resolved halo masses, i.e. for more than two decades in 
mass, in our simulation without outflows. This is the case 
even though a large fraction of the infalling gas in halos with 
Mh > 10^^'^ M0 is hot, reflecting the minor contribution of 



hot mode accretion to star formation at this redshift ( Keres 
|et al.|2009a] . 

In our constant-velocity wind simulations, 
SFR/ {Mcoid + Mism) is heavily suppressed in low- 
mass halos at all redshifts, as outflows eject gas from 
galaxies before it can form stars. In sufficiently massive 
halos, the situation is reversed and the star formation rates 
in halos can exceed cold gas-l-ISM accretion rates by factors 
of a few (see also Fig. [T]|. In all cases, we find that the 
SFR/ (McoM + Mism) ratio is enhanced by a factor < (1-f 77), 
which is consistent with the additional star formation being 
fueled by the fall back of wind ejecta (in this regime. 



Fig. [4] shows that AfcoW + Mism is relatively unaffected 
by winds, so that variations in SFR/{Mcoid + Mism) 
reflect variations in SFR). As expected, the high- velocity 
winds (fwinds) contribute enhanced star formation only 
in halos significantly more massive than the slow winds 
(winds, swinds) do. At late times {z = 0), the ratio 
SFR/{McoM + Mism) is suppressed much more heavily 
than the factor l/{l + r]) that would result if the only effect 
of outflows were to remove rj Mq of galactic gas for reach 1 
M0 of stars formed. This again indicates that outflows are 
able to prevent accretion of gas onto galaxies in this regime. 
Interestingly, the predicted z = star formation rates are 
lower than observationally-inferred ones in all of o ur wind 
models for Mh < lO^'^ Mq. For example, at z — 0.1, Conroy 



in Mh = 10^* M0 halos, and ~ 0.2 Mq yr" 
M0 halos, based on abundance matching. At z 
simulations predict comparable star formation rates in 
Mh ~ 10^^ Mq halos for all but the most extreme fwinds 
outflow implementation. For Mh « 10^^ Mq halos, however, 
the winds, swinds, and fwinds simulations predict only 
SFR 0.04, 0.02, 0.002 M© yr~\ respectively. 

The dominant factor explaining this discrepancy is 
likely the suppression of superwinds from real galaxies at 
z ~ 0. Observations show that powerful superwinds are 
ubiquitous in galaxies with star for mation rates per unit 



& Wechsler (|2009|) infer star formation rates ~ 3 M© yr 

Mn Mh = 10" 
0, our 



area Sspr ^0.1 Mq yr kpc (e.g., Heckman 20021), with 
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Figure 8. Comparison of the cold and hot halo gas fractions, for the different wind prescriptions. Solid curves: no winds. Dashed: 
constant-velocity winds with Dw = 342 km and mass loading rj = 1 (winds). Dotted: constant-velocity winds with Vv, = 342 km 
and mass loading r] = 2 (swinds). Dash-dotted: constant- velocity winds with Dw = 684 km and mass loading rj = 2 (fwinds). Vertical 
markers indicate the halo mass at which the cold and hot gas mass fractions are equal for the no-wind case in each panel. 



outflows from other galaxies being significantly weaker (e.g., 
Chen et al.|2010 1. Theoretically, this critical surface density 
can be derived from the condition that the massive star clus- 
ters in which most of the star formation of a galaxy occurs 
must be able to accelerate gas to velocities exceeding the 
escape velocity of the galaxy as a whole in order to produce 
a galactic-scale outflow ( [Murray et al.p011[ ). Whereas this 
condition is satisfied in local starbursts and high-redshift 
Lyman break galaxies (e.g., Erb et al. 20061, it is usually 



our simulations aside for instantaneous gas recycling from 



not in present-day ordinary spirals, like the Milky Way. The 
high star formation rates permitted by the large cosmolog- 
ical infall rates at high redshift (Fig. [l| ensure that super- 
winds operate effectively at early times. As the universe ex- 
pands and becomes more tenuous, however, ordinary galax- 
ies can no longer maintain star formation rates sufficient 
to drive significant outflows, fn our simulations, this condi- 
tion is not implemented and all galaxies drive outflows. As a 
result, the star formation rates at z ~ are artiflcially sup- 
pressed by feedback in the simulation with winds in some 
regimes. Note that star formation rates of the order of the 
observed ones are predicted in our simulation without winds 
at Mh ~ 10^^ — fO^^ Mq, indicating that cosmological infall 
is indeed sufficient to account for most of that star formation 
in the absence of feedback. Stellar mass loss, especially from 
AGB stars which can return up to ~ 50% of the mass of a 
stellar population into the ISM over long time scales, is also 
likely to permit modestly higher star formation rates at late 
time (e.g., Leitner fc Kravtsov|201l" l, but is not modeled in 



Type II supernovae at the 10% level (p^. The high ISM 
metallicities that would result at z = if most star forma- 
tion were fueled by stellar mass loss are however in tension 
with the measured ones, which argue for substantial recent 



accretion of metal-poor intergalactic gas (e.g., Sancisi et al 
|2008j . 

In high-mass halos, the situation is reversed and our 
simulations overpredict the observed star formation rates. 
This is the same problem alluded to in §3.11 with respect to 
the baryon mass function, and is indicative of the need for 
additional feedback, most likely from AGN. 

Another interesting feature of Figure [S] is that the ra- 
tio SFR/ {Mcoid + A^ism) is significantly below unity even 
without feedback in halos of Afh < lO'^^'^ M0 at z = 5, even 
though these halos are well into the cold mode-dominated 
regime. The likely explanation at this redshift, when the 
Universe just over 1 Gyr old, is that the Universe is so young 
that the recently formed low-mass galaxies have not had 
time to settle in a steady state in which the star formation 
rate follows the supply rate. In fact, the molecular gas de- 
pletion time scales, while yet to be directly measured at this 
redshift, are comparable to this time scale at 2 ~ 2 (e.g., 
[Genzel et aT||2010[ ). As Figure [5] indicates, an approximate 
steady state between cold gas supply and star formation rate 
is reached by z = 3 in our simulation without outflows. 

At 2 = 0, the star formation rate is again slightly sup- 
pressed relative to the influx of cold gas through the virial ra- 
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dius. Multiple effects likely contribute to this. First, present- 
day low-mass galaxies are particularly prone to support ex- 
tended gaseous discs, in which material accreted from the 
IGM can be temporarily stored, rather than rapidly con- 
verted into stars. At the low densities of outer discs, star 
formation is observed to become very i nefficient, with de - 
pletion time scales up to ~ 10^^ yr (e.g., Bigiel et al.|20ld |. 
The SFR therefore naturally deviates from the supply of cold 
gas when a substantial fraction of the infalling gas joins the 
outer parts. At larger masses, deviations between the cold 
gas-|-ISM accretion rate at the virial radius and the star for- 
mation rate within the halo are not surprising at z — 0. 

To this discussion we must add that the star formation 
rates in our simulations depend to some extent on the sub- 
resolution star formation prescription {[2.1). Furthermore, 



they are sensitive to the gas distribution in galaxies, such as 
the fraction of gas stored in low-density outer parts (Keres et 
al., in prep.). The quantitative details of the star formation 
rate predictions should therefore be treated with appropriate 
caution. The general points that the star formation rates in 
galaxies are typically supply-limited and fueled by cold gas 
accretion, and the importance of the re-accretion of wind 
ejecta, are however robust. 

3.2.3 When Does SFR ^ {nb/QoM)MDM ? 

The (Mcoid + Mism) - Mdm and SFR - (Mcoid + Mism) 
relations discussed above combine to yield the SFR — Mdm 
relation. This relation most directly quantifies the validity of 
the hypothesis that galaxies form stars at rates determined 
by the cosmological infall predicted in dark matter models. 
As Figure ^{SFR/MoM vs. halo mass and redshift) shows, 
the assumption that SFR « (r2b/f2DM)MDM is essentially 
never very accurate, even in simulations without galactic 
outflows. When outflows are included, the relation is sub- 
stantially modified. In the constant-velocity models explored 
here, the two principal effects are to suppress the star for- 
mation rates in low-mass halos, but to boost them in more 
massive ones that can confine their outflows and later fuel 
their galaxies via re-accretion. The resulting SFR — A/dm 
relation is however very sensitive to the uncertain outflow 
parameters, thus limiting the accuracy of simple estimates 
based on the dark matter infall. 



3.3 Integrated Baryonic Halo Properties 

We now turn to quantifying some of the integrated bary- 
onic properties of halos. Since all the particles within the 
halos are included in these predictions, we relax the mini- 
mum number of dark matter particles to 64 (from 500 for 
the accretion rates) to cover a larger dynamic range. 

Figure [7] shows how the baryonic mass in halos is di- 
vided between cold gas, hot gas, stars, and ISM, for our 
different simulations. In the no-wind case, the total baryon 
fraction (shown by the solid green curves) is close to the uni- 
versal ratio Qh/^m above a redshift-dependent halo mass, 
corresponding to suppression by the UV background. The 
order of magnitude of this mass scale corresponds to the 
IGM temperature Tigm ~ 10** K, multiplied by a factor of a 
few to account for adiabatic heating of the gas as it falls into 
halos, and is most directly related to the halo circular veloc- 
ity (e.g., |Efstathiou||1992[ [Thoul fc Weinberg||1996[ ). |Thoul 



& Weinberg ( 1996 1 showed that halos with circular veloc- 
ity ^ 50 km s~^ can maintain pressures high enough 
to prevent the gas from collapsing with the dark matter. 
Defining Vc = {GM^/ R^i.Y''^ and adopting Mh = Migot, 
7?vir = Ri&ob for consistency with the rest of this paper 



(5 2.3 1 



Mh = 10^° Mc 



V50 km s-i/ 



1 + g 
4 



-3/2 



(11) 



This result shows how the UV background suppression mass 
scale should increase with time, as is indeed found to be the 
case in our simulations. We note that this crude analytic 
derivation is only meant to illustrate the principal effects 
seen in Figure [7] In reality, the details of this suppression 
are also sensitive to the spectrum of the UV background 
via its impact on the photoheating rate and on the cooling 

and is subject 



2009 1 



function (e.g., Faucher-Giguere et al.| 
to a time delay ( |Gnedin|2000| . 

In the no-wind case, the total baryon fraction actually 
exceeds the universal value ilb/^^m by a small amount in 
some regimes. In |3.2.1[ we found a similar effect, but for 
the accretion rates. It is now apparent in Figure [7| that the 
super-universal baryon fractions occur in halos with mass 
just above the photoionization suppression mass scale, lend- 
ing credence to the hypothesis that the gas kept out of 
low-mass halos by photoheating provides the material neces- 



sary for the super-universal baryonic accretion rates (Keres 
|2007[ ). In very massive halos, gas pressure in hot atmospheres 
may push out a small fraction of the baryons outside the 
virial radius and explain the slightly sub-universal baryon 
fractions ( |Dave||2009| ). 

Figure [7] groups stars and ISM together to limit the 
number of curves shown. We conventionally refer to the com- 
bination as the "galactic mass", cautioning that the ISM 
boundary here only reflects the implementation of the mul- 
tiphase sub-resolution model (nn 5S 0.13 cm~^; see \ 2.A), 



and that gas stored in extended low-density discs may be 
missed by this measure. The relative contribution of the 
stellar and ISM components varies with both redshift and 
halo mass. At z = 0, > 90% of the galactic mass is in stars 
for all halo masses. This fraction gradually diminishes with 
increasing redshift, and preferentially so in low-mass halos, 
whose galactic mass is increasingly dominated by ISM gas. 
At z = 3, the stellar and ISM components divide the galactic 
mass equally at Afh « 10^'^'^ M©, with more massive halos 
having a greater mass in stars, and lower-mass halos having 
a greater mass in ISM. By z = 5, the 50% point is 0.5 dex 
higher, at Mh ~ 10^°'^ M©. Although these numbers are for 
the no-wind simulation, the picture is qualitatively similar 
in the simulations with outflows. 

As expected from their effects on the accretion rates, 
the principal impact of outflows is to remove baryons from 
the halos, especially at low masses. At very early times, the 
impact of winds in low-mass halos is small in our simula- 
tions (e.g., at Mh < 10^° Mq in the z = 5 panel of Fig. [7| 
because the young galaxies they host have not had time to 
form a lot of stars and to eject much material. By z = 0, 
the suppression of the baryon fraction in low-mass halos is 
however dramatic. At sufficiently high masses, halos con- 
fine the winds and so the predicted total baryon fractions 
approach Slb/f^m, as in the no- wind case. In general, the 
baryon fraction is increasingly suppressed as either the mass 
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loading factor 77 or the wind velocity «„ is increased. At fixed 
halo mass and outflow parameters, winds are also most ef- 
fective at suppressing the halo baryon fractions at low red- 
shift. This occurs for the same basic reason explaining the 
enhanced efficiency of preventive feedback from photoion- 
ization. Namely, halos of a given mass are less dense at low 
redshift, which reduces both their escape velocity and the 
hydrodynamic drag they can provide. 

In terms of observations, the compilation of |Dai et al.| 
| |2010[ ) suggests that the actual halo baryon mass fractions 
are close to universal only in clusters with mass Mh > 10^** 
^0 [i^c 600 km s~^) locally, and are increasingly sup- 
pressed below this mass. For halos below approximately this 
mass, which do not emit significantly in the X-rays, it is how- 
ever generally not possible to directly measure the diffuse 
halo gas. The estimated baryon mass fractions are therefore 
principally based on galactic baryonic mass measurements, 
and part of the baryons may be "missing" only because of 
observational limitations. For most galaxy-sized halos, how- 
ever, the detected baryons represent such a small fraction 
of the dark matter mass (e.g., McGaugh 2005 Hoekstra 
|et al.|2005| , and the halo cooling times should be sufficiently 
short, that it is unlikely that all the missing baryons are sim- 
ply undetected. In contrast, in our no-wind simulation, all 
halos with Afh ^ 10^^ M© have near-universal baryon frac- 
tions. Only our most extreme f winds simulation comes close 
to suppressing the baryon fraction significantly in halos as 
massive as Mh ~ 10^'' M© at z — 0, again highlighting the 
need for very strong feedback. This simulation is however 



too effective at removing baryons from L* galaxies (S2.2I, 
indicating that realistic feedback must have a more complex 
mass dependence, and is likely to have multiple sources. |Mc-| 
Carthy et al. (20101, for example, showed that the inclusion 
of AGN feedback, in addition to stellar feedback, can ac- 
count for the suppression of the baryon mass fraction in 
groups with mass up to Mh ~ lO^'' M0 . 

Figure[8]considers only the diffuse halo gas (i.e., exclude 
ISM material) and compares the mass fractions of the cold 
and hot components, as a function of halo mass and redshift. 
In agreement with previous work, the diffuse gas mass in 
low-mass halos is completely cold-dominated, and it is com- 



pletely hot-dominated in high-mass halos (e.g., Keres et al. 
|2005[[Keres et al.|2 009a). In the no- wind case, the halo mass 
at which the cold and hot halo gas mass fractions are equal, 
the transition mass Mt, agrees well with the equivalent but 
lower-resolution runs of [Keres et al.| ( |2009a[ ) at z = and 
2 = 3, with a value slightly below Mh = lO"'^ Mq. The 
z — 4: — 5 panels clearly show that Mt increases to higher 
values at higher redshifts, as the cooling times (particularly 
along the dense dark matter filaments feeding rare halos) 
are sufficiently short to permit the existence of cold streams 



even in halos that do develop virial shocks (e.g., Keres et al 
2005i [Keres et al.||2009a[ [Dekel fc Birnboim[[2006[ [Ocvirk 



et al.|2008 1 . A novel aspect of the present work is an explo- 
ration of several outflow parameterizations (see also [van de [ 
Voort et al.[[2011| ). At 



z = 0, the cold gas mass fraction in 
halos is seen to be nearly independent of the wind model, 
but some differences are seen at intermediate redshifts, par- 
ticularly near z = 2, where Mt in our most extreme outflow 
model (fwinds) is higher by more than 0.5 dex relative to 
the no- wind case. The main effect of outflows in our simula- 
tions here seems to be to increase the cold gas mass fraction 
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Figure 9. Net cold gas-|-ISM accretion rates through the virial 
shell at z = 0, 1, 2, 3, 4, 5, and 6 (from the bottom up). The solid 
grey lines show a fit to the 2^2 data, where the star formation 
rates are most directly connected to the cold gas infall rates. 



in halos by injecting them with cold gas from the ISM, and 
to consequently reduce the hot gas mass fraction. 

Although the transition mass in terms of mass fractions 
is approximately constant. Figure [l] shows that the equiv- 
alent transition point for the accretion rates, i.e. the mass 
where Mcom = Mhot at the virial radius, does evolve signif- 
icantly with redshift (see also Ocvirk et al.| 2008 ) . For some 
applications, such as in semi-analytic models which require 
the fraction of the gas accreted by halos in different forms, 
this redshift-dependent transition point may be more rele- 
vant. 



3.4 Fitting Formulae 

While the results of this paper highlight the sensitivity of the 
baryonic accretion rates to the feedback model, the accretion 
rates in the no-wind case do not depend directly on model 
assumptions and may therefore be generally useful, provided 
that the effects of feedback are separately accounted for. Of 
these, the cold gas accretion rate is the most important as it 
provides the main fuel for star formation. As shown in |3.1[ 
the infalling cold gas at the virial radius of massive halos is 
significantly suppressed during the infall aX z <2. The star 
formation rates also differ significantly from the cold gas 
accretion rates at the virial radius at 2 ~ 0. We therefore fit 
a simple model to the net AfcoW + Mism data versus halo 
mass and redshift, but limit ourselves to 2 2, where it may 
be most robustly used to model star formation. 

As it is useful to also consider the efficiency with which 
dark matter channels cold gas into halos, we begin by pro- 
viding a fit to our median dark matter accretion rates (Fig. 
[T|. This fit accounts for the minor differences with the corre- 
sponding FMBIO formula arising from the definition of our 
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accretion rates as being instantaneous and for the slightly 



different cosmology (^3.1 1. Whereas the Millennium simula- 



10" 



tions analyzed by FMBIO contained only dark matter, ours 
explicitly include baryons. Our dark matter accretion rates 
therefore, by definition, differ by a factor {Q,ni — r2b)/fim- 

We adopt the same functional form that was shown to 
yield an accurate fit by FMBfO, 

{A/)^ed = N{1 + az) " V^m{l + zr + QA. 

(12) 

For the dark matter only, we find best-fit values (A'^, a, /3) = 
(33.6 M0 yr"\ 0.91, 1.06) to the data at z = 0, 1, 2, 3, 4, 5 
and 6 simultaneously. These fit the data in Figure [l] nearly 
perfectly. For the cold gas+ISM data, we fix the redshift 
evolution parameter a to the best-fit dark matter value, but 
allow for an additional overall redshift tilt factor (1 + z)'' . 
The resulting fit is as good as leaving a free, but yields 
a cleaner result for the efficiency below. The best-fit pa- 
rameters for (Mcoid + AfisM)mcd {z ^ 2 Only) are then 
(iV, a, /?, 7) = (1.84 Mq yr~\ 0.91 0.81, 0.38). 
These fits allow us to evaluate the efficiency 



Ecold 



mcd 



{Mum ) mcd [^^b / ( n 



(13) 



which would be unity if halos were supplied with cold 
gas-|-ISM with a universal baryon fraction of the dark matter 
accretion rate. Using the above fits. 



0.47 



1 



1012 Mp 



(14) 



Overall, the cold gas-|-ISM accretion rate has a shallower de- 
pendence on halo mass than the dark matter accretion rate. 
The owes principally to the greater ability of massive halos 
to main stable virial shocks, which heat a fraction of the gas 



(S3.3I. These new fits to the cold gas-|-ISM accretion rates, 
which primarily drive star formation, provide a more real- 
istic basis for galaxy formation models, and may have im- 
plications for studies that have assumed steeper mass slopes 
based on dark matter results (e.g., Bouche et al.|2010 l. Note 
that equation ( 14 1 would (unphysically) predict ecoid — > oo 
as z — >■ oo or A/h — >■ 0, so that the fits presented here should 
be understood to be valid only over the parameter space 
covered by the numerical data shown in Figure |9] at the 
precision that the fit to (Afcoid + A^isM)mcd represent those 
data. 

Finally, it is instructive to compare Figure |9] (net 
Afcoid + AfisM) to Figure [2](infaning 

Afcoid component only). 
The similarity of the curves at corresponding redshifts indi- 
cates that the net cold gas accretion rates, in the no-wind 
simulation, are approximately equal to the infalling com- 
ponent only. Furthermore, the ISM component is a sub- 
dominant component of A/cow + A/ism. 



4 DISCUSSION AND CONCLUSION 

In this paper, we systematically quantified the rates at which 
dark matter halos accrete baryons and the resulting baryon 
mass fractions using a set of four high-resolution cosmolog- 
ical hydrodynamical simulations. Motivated by recent work 
indicating that star formation in galaxies is primarily driven 
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Figure 10. Comoving number density of halos as a function of 
minimum net cold gas-|-ISM accretion rate through the virial shell 
in our no-wind simulation at 2 = 0, 1, 2, and 4 (z = 3 omitted 
for clarity) . The curves terminate at the highest accretion rates in 
the 40 comoving Mpc box. For reference, the star shows the 
accretion rate equal to the median star formation rate in the Ha 
sample of |Erb et al. | | |2006[ l and the estimated the number density 
of BX galaxies based on clustering ( [Adelberger et al.|2005| . 



by the supply of cold gas, we computed the gas accretion 
rates for the diffuse cold gas (nn < 0.13 cm~^, T < 2.5 x 10^ 
K), hot gas (riH < 0.13 cm-^ T > 2.5 x lO'^ K), and ISM 
(tih 0.13 cm~"^) components separately, in addition to 
the dark matter. A main goal was to understand how ex- 
isting results on the growth of the dark matter component 
of halos, which is accurately described by a set of analytic 
tools, including the Press-Schechter and excursion set for- 
malisms ( Press & Schechter 1974 Bond et al. 1991 1 and 



fitting formulae for their merger and accretion rates (e.g., 
Wechsl er et al.|[2002l [Neistein fc Dekcl, 2008, Fakhouri S 
Ma||20 08 ; Fak houri et al.||2010| [McBride et al.p OOO 



Genel 



et al. 2010|, can be extended for baryons. Currently, dark 



matter-based analjd;ic tools form the basis of many theoret- 
ical studies of galaxy formation, but it is usually necessary 
to make important assumptions to apply them to baryons. 

We investigated the baryonic accretion rates as a func- 
tion of halo mass and redshift in a simulation without strong 
feedback, and for three otherwise identical simulations that 
implemented a phenomenological outflow model with dif- 
ferent parameters (Fig. [T|. This allowed us to explore the 
sensitivity of our results to the uncertain feedback physics. 
By focusing on the infalling component only, we showed that 
winds can not only suppress the net accretion rates of cold 
gas by transporting cold gas outward, but also substantially 
affect the infalling cold material hydrodynamically, partic- 
ularly in low-mass halos at low redshift (Fig. |2|. In halos 
sufficiently massive to confine their outflows owing to a com- 
bination of gravitational deceleration and hydrodynamic in- 
teractions, the accretion rates through the virial radius are 
however robust to the wind prescription. By evaluating the 
cold gas accretion rates through radial shells of radii ranging 
from 1-Rvir to 0.2i?vir (Fig.|3|, we explicitly showed how the 
accretion rate of cold gas diminishes as it penetrates into ha- 
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los, owing to a combination of shocking and accretion onto 
satellite galaxies. Our results also confirm the persistence 
of cold streams into massive halos with hot atmospheres 
at z ^ 2 (IKeres et al.||2005| |Keres et al.||2009a| |Dekel & 



Birnboim|2006 Ocvirk et al.|200"8 1 and their gradual disap- 



pearance a lower redshift. 

We directly tested the validity of a number of common 
simplifying assumptions, including whether the cold gas ac- 
cretion rate traces the dark matter accretion rate with a uni- 
versal baryon mass fraction (Fig. [4|, and whether the star 
formation rate within halos follows the cold gas accretion 
rate at the virial radius (Fig. [5|. Ultimately, we quantified 
the accuracy with which the star formation rates in halos 
can be predicted from their dark matter accretion rates (Fig. 
[6|. We find that in general such assumptions are not very 
accurate, highlighting the need for more accurate hydrody- 
namical results. In the absence of galactic outflows, the dark 
matter accretion rate is a predictor of the star formation rate 
at the factor of ~ 2 level in halos below the transition mass 
Mt ~ 10^^'^ Mq between the cold and hot modes. Above 
this transition mass, the dark matter accretion rate is an 
increasingly poor predictor of star formation, as the cooling 
times become longer than the free fall times, and gas enter- 
ing the virial radius is no longer efficiently channeled into 
galaxies. This is especially so at lower redshift. 

When galactic outflows are included, the star formation 
rates are sensitive to the outflow parameters, especially in 
low-mass halos from which gas can be ejected altogether, 
where in some regimes star formation can be effectively 
stopped. In massive halos confining their outflows, the in- 
stantaneous star formation rates can on the other hand be 
boosted relative to the no-wind case by the re-accretion of 
gas onto galaxies. 

For our no-wind simulation, which is relatively free of 
ad hoc assumptions, we have provided simple fitting formu- 
lae for the infall rates of cold gas through the virial radius, 
and for the corresponding efficiency with which dark matter 
channels cold gas into the halos (eq. 1 12]|l4 l). These will per- 
mit improvements of existing studies based on dark matter 
only. 

An interesting consistency check is to compare the mea- 
sured star formation rates at high redshift, where the high 
values are believed to be allowed by enhanced intergalactic 



accretion rates (e.g., Dekel et al. 2009 1, to the accretion rates 



of cold gas predicted by our simulations. We focus specifi- 
cally on the 2 ~ 2 star-forming galaxies that have been the 



subject of intense observational studies recently (e.g., Stei- 
del et al.||2004| |Erb et aH20()6| |Daddi et al.|2007| ), and the 



targets of detailed integral field spectroscopic (IFS) analyses 



(e.g.,lGenzel et al.|2006 


Forster Schreiber et al.||2009 Law 


et al.|2009 Wright et al. 


20091. These IFS studies have re- 



vealed the 2; ~ 2 disc galaxies to be very clumpy, which has 
been argued to be a direct consequence of high gas fractions 



maintained by rapid accretion (e.g., Elmegreen & Elmegreen 
[2006} [Dekel et al.|[2009{ [Genzel et al.||2011[ ). An exact com 
parison is complicated by the complex and sometimes un- 
even observational selection, and because the distribution 
of halo masses hosting the galaxies is not fully captured by 
the characteristic halo mass inferred from clustering. Fur- 
thermore, significant uncertainties exist in star formation 
rate estimates. For example, recent calorimetric SFR mea- 
surements by the Herschel Space Observatory suggest that 



previous UV-based SFR estimates were biased high by a fac- 
tor ~ 2 (Nordon et al. 20101. In general, satellite galaxies 



can also accrete material infalling through the virial radius 
of main halos, so that the halo star formation rates will be 
shared between the central galaxies and satellites. Keeping 
these caveats in mind, we therefore consider only the broad 
picture here. 

The median star formatio n rate in the IFS sample of 
Forster Schreiber et al] (|2009| is SFR « 70 Mq yr"^ The 



clustering of galaxies in a similarly selected parent sample 
indicates a characteristic halo mass Mh = 10^^'* — lO^'^ '^ 
M0 (Adelberger et al. 20051. On the other hand, our fit 



to the z 2 simulation data ( ^3.4[ ) implies that the cold 
gas accretion rate through the virial radius, including ISM 
material, is ^ 23 M© yr"^ for a Mh = 10^^ Mq halo at 
z — 2 in the absence of outfiows. While this is short by 
a factor ~ 3 to explain the median SFR ~ 70 Mq yr~^, 
the selection criteria of the IFS sample favor above-average 
systems permitting sp atially-resolved spectroscopy ( |Forster| 



Schreiber et al. 



20091 



The 90'" percentile Mcow + Mism 
accretion rate for this halo mass and redshift is in fact a 
factor ~ 3 higher than the median. 

An alternative comparison is with the Ha sample of 



Erb et al. (20061, which was d rawn from the large r rest 
UV-selected sample described by Steidel et al. ( 2004| p] and 
which has a comparable characteristic halo mass Mh ~ 10^^ 
Mq. The median SFR Ri 30 Mq yr"^ of this sample is in 
good agreement with the predicted supply rate. A potential 
worry is that the outflows inferred in these systems, with 
mass loading factors up to a few times their star formation 
rates (e.g., Erb et al. 2006| [Erb|pC)08j , might require com- 
mensurately higher accretion rates to sustain the star for- 
mation. However, as Figures[T]and|6]show, the halos are suf- 
ficiently massive that, for realistic parameters, outfiows can 
actually boost the star formation rates via gas re-accretion. 
Fig ure |10| also shows that the number density of BX galax- 
ies ( Adelberger et aL]|2005 l, from which Erb et al.'s sample 
is drawn, agrees well with the number density of halos with 
cold gas-flSM accretion rates sufficient to fuel their star for- 
mation in our no-wind simulation at z = 2. Gas recycling 
from supernovae and AGE stars can also sustain star for- 
mation rates exceeding the external supply rate by ~ 50%, 
and hot mode accretion can provide some extra fuel (e.g., 
Keres et al.pOOS , Kcrcs et al 2009aj) 



While these considerations show that the high SFR ~ 
10 — 100 Mq yr~^ of z ~ 2 galaxies are, broadly, well ex- 
plained by being primarily determined by cold mode ac- 
cretion, the quantitative details should be interpreted with 
caution. In particular, typical samples (e.g., Erb et al.|2006 



Daddi et al.]|2007 l cover up to two orders of magnitude in 



^ UV selection techniques are designed to discover star-forming 
galaxies and therefore miss passive galaxies (e.g., [van Dokkum et| 
|al.|2006"| . This population may account for as much as half of the 
massive galaxies at 2; ~ 2, but their very low star formation rates 
(e.g., |Kriek et al.|2006[ l suggest that processes not included in our 
simulations (such as feedback during galaxy mergers) terminated 
star formation in them. Provided that, at a given halo mass, this 
putative quenching process is stochastic in the sense that it does 
not depend significantly on the cold gas accretion rate, it is fair 
to compare our simulation results with the sub-sample of star 
formers. 
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stellar mass, so that comparisons with characteristic halo 
masses are too simplistic to draw firm conclusions. The more 
intense star- formers with SFR > 100 Mq yr~^ (including 
sub-millimeter galaxies; e.g., Hayward et al. 2011) are most 
likely possible owing to a combination of variance in the 
accretion rates and other fueling channels, such as gas re- 
cycling, hot mode accretion, and mergers|^A more detailed 
comparison of these with the results of simulations would 
be an interesting topic for future work, but is challenging 
owing to observational uncertainties. 

At 2 = 0, we have shown that cosmological accretion 
can also account for the star f ormation rates inferred in 
Mh ~ 10" -10^2 Mq halos (e.g., |Conroy fc Wechsler|2009 1. 
Our simulations with outfiows however all under-predict the 
star formations rates at Mh ~ lO'^^ M©, by factors ranging 
from 5 to 100. Most likely, this is because whereas all star- 
forming galaxies at high redshift drive powerful superwinds, 
this is not true at z ~ 0. Instead, it appears that only galax- 
ies with star formation rate surfa ce density Ssfr 55 0.1 M0 
yr~^ kpc~'^ have superwinds (e.g., Heckman|2002 I. This cri- 
terion is usually satisfied at high redshift, where the star for- 
mation rates are high, but locally only in starbursts. Such a 
critical star formation rate surface density may be related to 
the ability of feedback processes to eject gas from the host 
galaxy (e.g., Murray et al.|20TT l and is not modeled in our 
phenomenological outflow implementation. It is also miss- 
ing in parametric momentum-driven wind scalings that do 



sition point in the accretion rates at the virial radius does 
depend significantly on redshift. 

Finally, these results should be useful to guide the im- 
plementation of cold accretion in semi-analytic models of 



not self-consistently model the wind generation (e.g.. Op- 



penheimer & Dave 2008 Oppenheimer et al. 20101. This 



suggests that more realistic wind models will be necessary 
to reproduce the observed local galaxy population. 

Without winds, the total baryon fraction in halos is 
close to universal in all halos above the UV background pho- 
toheating suppression scale, which our simulations resolve at 
all redshifts z < 5 (Fig. [7| . The UV background suppression 
mass scale is seen to shift to progressively higher masses with 
time (from ~ 10^° Mq at ^ = 3 to ~ 10^"'^ Mq at z = 0), 
refiecting the fact that the most relevant halo property is 
the circular velocity, which increases with redshift at fixed 
halo mass. Such universal baryon fractions are in contrast 
with observations, which suggest that only galaxy clusters 
with Mh > 10^* M0 may have near-universal baryon frac- 
tions at z ~ (e.g., Dai et ar]|2010 but see the discussion 



of the uncertainties on this mass scale in S3.3I. Outfiows 



are effective at suppressing the total baryon fractions up to 
a halo mass depending principally on the outfiow velocity. 
Only our most extreme fwinds outfiow prescription, with 
Vv, = 684 km s~^, suppresses the baryon fraction in halos 
close to the necessary scale. This again underscores the need 
for strong feedback to match the observations, likely from 
AGN in addition to stars. 

While the transition halo mass Pit ~ 10^^'^ Mq at 
which halos host an equal mass fraction of cold and hot 
gas is approximately constant with redshift and insensitive 
to the outfiow prescription (Fig.[8|, the corresponding tran- 



^ During major mergers, gas stored in the galaxies can be effi- 
ciently consumed in a nuclear starburst, resulting in an instanta- 
neous star formation rate tliat can greatly exceed the cosmological 
infall rate (e.g., Mihos &: Hernquist|1996 l. 



galaxy formation (e.g., Cattaneo et al. 2006 Khochfar & Silk 
| 2009[|Lu et al.||2011[ [Benson fc Bower||2011[ ), in particular 
in understanding how galactic outfiows modify the standard 
results from hydrodynamical simulations. So far, the wind 
models we have explored have been phenomenological and 
the parameters were chosen to cover a representative range 
of possibilities. In future work, we will study the properties 
of self-consistently generated outfiows ( [Hopkins et al.|2011[ ), 
and the details of their interactions with infalling material. 
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